Dynamics of the Magnetic Flux Trapped in Fractal Clusters 
of a Normal Phase in Percolative Superconductors 
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The effect of the fractal clusters of a normal phase, which act as pinning centers, on the dy- 
namics of magnetic flux in percolative type-II superconductor is considered. The main features 
of these clusters are studied in detail: the cluster statistics is analyzed; the fractal dimension of 
their boundary is estimated; the distribution of critical currents is obtained, and its peculiarities 
are explored. It is found that there is the range of fractal dimension where this distribution has 
anomalous statistical properties, specifically, its dispersion becomes infinite. It is examined how the 
finite resolution capacity of the cluster geometric size measurement affects the estimated value of 
fractal dimension. The effect of fractal properties of the normal phase clusters on the electric field 
arisen from magnetic flux motion is investigated for the cluster area distribution of different kinds. 
The voltage-current characteristics of fractal superconducting structures in the resistive state are 
obtained for an arbitrary fractal dimension. It is revealed that the fractality of the boundaries of the 
normal phase clusters intensifies the magnetic flux trapping and thereby raises the critical current 
of a superconductor. 

PACS numbers: 74.81.-g; 74.25. Qt; 74.25. Sv, 74.25.Fy; 74.81.Bd 



I. INTRODUCTION 

An essential feature of the clusters of columnar defects 
in superconductors consists in their capability to trap a 
magnetic fluxpiiSi 3 - Holding in place the vortices driven 
by the Lorentz force, such clusters can act as effective 
pinning centers. 4,5 This feature can be used in making 
new composite superconducting materials of enhanced 
critical current valuei&L^ The characteristics of the nor- 
mal phase clusters exert an appreciable effect on vortex 
dynamics in superconductors, especially when the clus- 
ters have fractal boundaries i 9 ' 10 ! 11 ' 12 ! 13 In the present pa- 
per the magnetic and transport properties of composite 
superconductors with fractal clusters of a normal phase 
will be reviewed as well as the phenomena limiting the 
current-carrying capability of such superconductors will 
be considered. 

Let us consider a superconductor containing columnar 
inclusions of a normal phase, which are out of contact 
one with another. These inclusions may be formed by the 
fragments of different chemical composition, as well as by 
the domains of reduced superconducting order parame- 
ter. The similar columnar defects can readily be created 
during the film growth process at the sites of defects on 
the boundary with the substrate. When such a supercon- 
ducting structure is cooled below the critical temperature 
in the magnetic field oriented along the direction of the 
longest size of these inclusions, the magnetic flux will be 
frozen in the normal phase clusters. By the cluster we 
mean a set of the columnar defects, which are united by 
the common trapped flux and are surrounded by the su- 
perconducting phase. Even after the external field has 
been turned off, the flux trapped in these clusters is kept 
unchanged due to the currents that are steadily circu- 
lating around them through the superconducting loops. 



The distribution of the trapped magnetic flux resulting 
from such a magnetization in the field-cooling regime will 
be two-dimensional. So, instead of dealing with an ex- 
tended object, which indeed the normal phase cluster is, 
we will consider its cross-section by the plane carrying a 
transport current. We will examine the geometric prop- 
erties of the normal phase clusters in the planar section 
only, where the boundaries of the clusters are statistically 
self-similar. 

The paper is organized as follows: The setting of the 
problem is described as well as the relationship between 
magnetic flux trapping, pinning properties, and geomet- 
ric characteristics of the normal phase clusters in super- 
conductor is discussed in Sec. 2. A geometric probability 
analysis of the weak links distribution over the cluster 
perimeter in terms of path integrals is made in Sec. 3. 
The concept of fractals and the basic notions in this field 
are reviewed briefly in Sec. 4. The statistical properties 
of normal phase clusters in YBCO superconducting films 
are studied, and the fractal dimension of their boundaries 
is estimated in Sec. 5. The voltage-current characteris- 
tics of superconductors with fractal clusters of a normal 
phase are obtained, and the peculiarities of the resistive 
state of the fractal superconducting structures are con- 
sidered in Sec. 6. 



II. MAGNETIC FLUX TRAPPING IN 
PERCOLATIVE SUPERCONDUCTORS 

The situation is sketched out in Fig. [IJa). The com- 
posite superconductor containing clusters of a normal 
phase represents a percolation system, where both the 
electric percolation of the supercurrent and the perco- 
lation of a magnetic flux may happen. If the part of 



the film surface covered by the normal phase exceeds the 
percolation threshold (which is equal to 50 percents for 
two-dimensional percolation 1 ^), the magnetic flux could 
be free to move in the transversal direction. There is 
no pinning at all in this situation, so it is out of our in- 
terest. Let us suppose that the opposite case is realized 
when the relative portion of superconducting phase ex- 
ceeds the percolation threshold (see Fig.[T](b)), so there is 
a superconducting percolation cluster in the plane of the 
film where a transport current can flow. Such a struc- 
ture provides for effective pinning and thereby raises the 
critical current, because the magnetic flux is locked in 
finite clusters of a normal phase, and so the vortices 
cannot leave them without crossing the surrounding su- 
perconducting space. Let us denote the total magnetic 
flux trapped in the superconductor after the field-cooling 
magnetization by <f>. If the transport current is passed 
through the sample, the trapped flux remains unchanged 
as long as the vortices are still held in the normal phase 
clusters. The larger clusters have to provide for a weaker 
pinning than smaller ones, because the larger the clus- 
ter size, the more weak sites, through which the vor- 
tices can pass under the action of the Lorentz force, can 
be located along its perimeter. When the current is in- 
creased, the magnetic flux will break away first from the 
clusters of smaller pinning force, and therefore, of larger 
size. Magnetic flux trapped into a single cluster is pro- 
portional to its area A. Therefore, the decrease in the 
total trapped flux A$ is proportional to the number of 
the clusters of area larger than a given magnitude A. 
So it can be expressed with the cumulative probability 
function W (A) = Pr{VA,- < ^4}, which is equal to the 
probability to find the cluster of area Aj smaller than a 
preset value of A: 



— = 1-W(A) 



(1) 



The left hand side of this formula is equal to the rel- 
ative decrease in the total trapped flux caused by the 
transport current of the same amplitude as the depinning 
current of the cluster of area A, and the right hand side 
gives the probability to find the cluster of area greater 
than A in the whole population. 

When the magnetic flux starts to break away from the 
normal phase clusters, the vortices must cross the sur- 
rounding superconducting space, and they will first do 
that through the weak links, which connect the normal 
phase clusters between themselves (see Fig. [ljb) where 
weak links are shown by the curves connecting the shaded 
regions, which designate the normal phase clusters). In 
Fig. [He) the equivalent weak-link-network circuit repre- 
sentation of such a system is shown. The normal phase 
clusters are presented by the cells of the net interlaced 
by the superconducting loops. Each of the loops con- 
tains weak links (shown by the crosses), which join the 
adjacent cells and so enable the vortices to pass from one 
cluster to another. 

Weak links form readily in high-temperature supercon- 





FIG. 1: Schematic drawing of superconductor with a mag- 
netic flux trapped in the clusters of columnar defects (a), its 
percolation representation in the cross-section by the plane 
carrying a transport current (b), and the equivalent weak- 
link- network circuit (c) . The sections of normal phase clusters 
are shaded. In picture (b) weak links are shown by the curves 
joining together the normal phase clusters, and, in picture (c), 
by the crosses on the superconducting loops that surround the 
normal phase inclusions. 



ductors (HTS) characterized by an extremely short co- 
herence lengthj 15 i 16 i 17 Various structural defects, which 
would simply cause some additional scattering at long co- 
herence length, give rise to the weak links in HTS. On a 
mesoscopic scale twin boundaries are mainly responsible 
for weak link existence ! 18 ' 19 ! 20 ' 21 Twins form especially 
readily in YBCO superconductors inasmuch as their unit 
cell is only close to the orthorhombic one. The twins 
can be spaced up to several nanometers apart, so even 
single crystal may have the fine substructure caused by 
twins. Magnetic and transport properties of HTS depend 
strongly on the orientation of twin planes with respect 
to the applied magnetic fiek h 22 ' 23 The flux can easily 
move along the weak links formed by twins i 24 ' 25 ' 26 ' 27 At 
last, on a macroscopic scale there are manifold struc- 
tural defects which can form weak links: that may be 
grain or crystallite boundaries as well as barriers aris- 
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FIG. 2: A simple illustration how the magnetic flux exits 
from a normal phase cluster through the weak link. Here 
It ~ Ii+ Ii is the transport current; is the current related 
to the trapped magnetic flux "I>; Fl is the Lorentz force, (a) 
- magnetic flux remains to be locked in the cluster as long as 
I2 + i* < Ic, where I c is the critical current of the weak link; 
(b) - as soon as I2 + > I c the Lorentz force expels the flux 
through the right weak link which has become resistive. 



ing from the secondary degrading the non-stoichiometric 
crystal into the domains with a high and low content 
of oxygen i 17 ' 28 i 29 ' 30 i 31 Moreover, a magnetic field further 
reduces a coherence length, thus resulting in more easy 
weak links formation^ In conventional low-temperature 
superconductors the weak links can be formed due to the 
proximity effect in sites of minimum distance between 
the next normal phase clusters. 

As soon as the transport current is turned on, this one 
is added to all the persistent currents, which maintain the 
magnetic flux to be trapped. The situation is sketched 
in Fig. [2j Each of these currents (like i$) is circulat- 
ing through the superconducting loop around the nor- 
mal phase cluster wherein the corresponding portion of 
the magnetic flux is locked (Fig.^a)). The loop contains 
weak links that join the adjacent normal phase clusters 
transversely to the path of the current. As the trans- 
port current is increased, there will come a point when 
the overall current flowing through the weak link will 
exceed the critical value, so this link will turn into a re- 
sistive state. As this takes place, the space distribution 
of the currents throughout the superconducting cluster is 
changed in such a way that the resistive subcircuit will be 
shunted by the superconducting paths where weak links 
are not damaged yet. Magnetic field created by this re- 
distributed transport current acts via the Lorentz force 
on the current circulating around the normal phase clus- 



ter. As a result, the magnetic flux trapped therein will 
be forced out through the resistive weak link, which has 
become permeable to the vortices (Fig. EJb)). So the 
cluster will not contain the magnetic flux any more. 

Thus, whatever the microscopic nature of weak links 
may be, they form the channels for vortex transport. It 
appears that according to their configuration each nor- 
mal phase cluster has its own value of the critical current, 
which contributes to the overall statistical distribution. 
By the critical current of the cluster we mean the cur- 
rent of depinning, that is to say, such a current at which 
the magnetic flux ceases to be held inside the cluster of 
a normal phase. When a transport current I is grad- 
ually increased, the vortices will break away first from 
clusters of small pinning force (weaker than the Lorentz 
force created by the transport current), and therefore, of 
small critical current. Thus the decrease in the trapped 
magnetic flux A<!> is proportional to the number of all the 
normal phase clusters of critical currents less than a pre- 
set value. Therefore, the relative change in the trapped 
flux can be represented by the cumulative probability 
function F (I) = Pr {VI j < 1} for the distribution of the 
critical currents of clusters: 



It is obvious that the critical current distribution 
F = F (I) is related to the cluster area distribution 
W = W (A) , because the cluster of a larger size has 
more weak links over its boundary with the surround- 
ing superconducting space, and thus the smaller current 
of depinning. 



III. GEOMETRIC PROBABILITY ANALYSIS 
OF THE DISTRIBUTION OF WEAK LINKS 
OVER THE CLUSTER PERIMETER 

In order to find out the relationship between the distri- 
bution of the critical currents of the clusters and the dis- 
tribution of their areas, the distribution of entry points 
into weak links over the perimeter of a normal phase 
cluster should be analyzed. The problem of exit of a vor- 
tex from a normal cluster represents the two-dimensional 
analogue of a problem of a random walk particle reach- 
ing a border i 33 i 34 At the same time, unlike the classic 
problem on the distribution of the exit points, here the 
boundary of the area is not absorbing all over. There are 
only discrete absorption points, which are located just at 
the sites where weak links are going out on the cluster 
boundary. In other words, these points are the points 
of the entry of vortices into weak links, or simply, en- 
try points. For simplicity, we will assume that after the 
vortex reaches the entry point, it passes all the way be- 
tween two adjacent normal phase clusters without being 
trapped inside the weak link itself. Here the magnetic 
flux is transferred by Josephson vortices. The Joseph- 
son penetration depth is large enough in the considered 
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materials, so the size of the region, where the vortex is lo- 
calized, much exceeds the characteristic length of all pos- 
sible structural defects that can occur along the transport 
channel. Thus the probability that such a vortex, driven 
by the Lorentz force, will be trapped in passing through 
a weak link is very small. This assumption agrees well 
with the results of research on the magnetic flux mo- 
tion along weak links2i^^ including twins ^^1^^^ 
At the same time, it allows us to highlight the role played 
by the cluster boundary in the magnetic flux dynamics. 

Let us consider the distribution of entry points over 
perimeter of a normal phase cluster. Generally, this dis- 
tribution varies from one cluster to another, so that each 
normal phase cluster has the entry point distribution 
function ip (I) of its own, which belongs to some func- 
tion class il. Here I is the co-ordinate measured along 
the cluster perimeter, so I S (0, P). In this context the 
functions of class Q are random elements of the statis- 
tical distribution. The probability distribution of func- 
tions ip (I) over all the clusters can be characterized by 
the functional Pr {ip (I)}, which is equal to the probabil- 
ity of finding a given function ip (I). 

In the most general way the geometric probability anal- 
ysis of the entry point into weak link distribution can be 
carried out by means of path integral technique.— The 
probability that the function ip (I) may be found within 
the class 51 can be expressed by the path integral 



Pr{iP (0| 0} 



Dip (l)Pr{ip(l)} 



Therefore, the most probable function of entry point 
distribution is the mean over all functions of class f2 



<S> (I) = iP (I) = / DiP(l)iP(l)Pr{iP(l)} 



(3) 



(£2) 



The path integral Fourier transform on the probability 
functional Pr {ip (I)} represents the characteristic func- 
tional 



H[k(l)] 



J ViP (I) exp (if dlk (I) ip (/)) Pr {ip (/)} 

/ ViP(l) Pt{iP(1)} 
(n) 



(4) 

where k — k (I) are the functions of a reciprocal function 
set, and integration in the kernel exp (i § dlk (I) ip (I)) is 
carried out over the cluster perimeter. 

The characteristic functional is the path integration 
analog for the usual moment-generating function. The 
probability functional Pr {ip (I)} can be written as an in- 
verse path integral Fourier transform on the characteris- 
tic functional 



Pr {iP (I)} = J Dk {I) exp l-i j dlk {I) ip (I) H [k (I)] 

where the path integration is carried out on the reciprocal 
function space. 



In the simplest case, when all the clusters are of 
equal entry point distribution, which coincides with the 
most probable one of Eq. ^ , the probability functional 
Pr{i/>(0} is zer0 f° r an "0(0 that, differ from ^ (/), 
whereas Pr {fy (I)} = 1. As this takes place the char- 
acteristic functional of Eq. ([J| becomes 



H [k {I)} = exp (i j> dlk (I) * (I) 



(5) 



If all entry points had fixed co-ordinates lj instead of 
the random ones, their distribution would be ip (I) = 

/3 Y)j—i S (I — lj) , where N is the number of entry points 
along the cluster perimeter, 8 (I) is Dirac delta function. 
The constant (3 is being chosen to normalize the distri- 
bution function ip (I) to unity, so that j3N = 1. 

Now suppose that all the points of entries into weak 
links are randomly distributed with uniform probabil- 
ity over the cluster perimeter, so the probability to find 
any j-th point within some interval dlj is proportional to 
its length. In that case the characteristic functional of 
Eq. ([J]) takes the form 



H[k(l)} 



N I N 

§ J] dlj exp \i(3 £ § dl k (I) S (I - l, 

3 = 1 \ 3=1 

N 

§ n 

■pN j n dl i ex p f *0 Yl k j = 

^Ujdlje^ 



(6) 



Expanding the function e l ^ k ^ in a power series at 
TV >> 1, and taking into account the condition (3N = 1, 
we may write 



j dle lf3k ^ = exp (i^ jdlk (I) 
that, after substitution into Eq. ([6]), gives 

H[k{l)]=exp(i-^- j dlk(l)J (7) 



The found characteristic functional of Eq. ([7]) has the 
form of Eq. (0 for the function of the uniform distribu- 
tion of entry points 



*(0 



f 

p 



(8) 



This means that all the clusters have the same uniform 
distribution of the entry points of Eq. (J5|) , for which the 
probability of finding a weak link at any point of the 
perimeter is independent of its position. 

Let us suppose that concentration of entry points into 
weak links per unit perimeter length n — N/P is constant 
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for all clusters, and all the clusters are statistically self- 
similar. In this case the mean number of entry points N 
along the cluster perimeter is proportional to its length: 



N = <pn(l)dl = nP 



(9) 



Next step will consist in finding the relationship be- 
tween the size of a cluster and its critical current. The 
pinning force corresponds to such a current at which the 
vortices start to break away from the cluster. As the 
transport current is increasing, the Lorentz force, which 
expels the magnetic flux, increases as well. The vor- 
tices start to leave the normal phase cluster when the 
Lorentz force becomes greater than the pinning force. 
At the same time, growing in current will result in 
re-distribution of the magnetic flux, which will pene- 
trate deeper and deeper into a transition layer on that 
side of the surrounding superconducting space where the 
Lorentz force is directed (see Fig.^b)). In order to leave 
the normal phase cluster, vortices have to reach the en- 
try points into weak links. The exit of the magnetic flux 
can be considered as the result of random walks of vor- 
tices driven by the Lorentz force, which is pushing them 
into weak links. A similar approach has been success- 
fully applied in analyzing the magnetic flux penetration 
in SQUID arrays^ 5 - The mean number of the entry points 
N available on the cluster perimeter provides the prob- 
ability measure of the number of the random walk out- 
comes, which are favorable for the vortex to go out. In 
the case of the uniform entry point distribution, from 
Eq. ([5]) it follows that N (x P, so the perimeter length 
also represents the probability measure of the amount of 
favorable outcomes for vortex to leave the cluster. The 
more entry points into weak links are accessible for ran- 
dom walk vortices, the smaller is the Lorentz force re- 
quired to push the flux out. Hence, we may write the 
following relationship between the critical current of the 
cluster and its geometric size: 



II II II II 
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r 1 1 

i a = a — 

N P 



(10) 



This expression is valid for the simplest case of uniform 
distribution of entry points, which is assumed to be the 
same for all clusters. Such a simplification allows us to 
emphasize that in the case being considered the magnetic 
flux is held in the normal phase cluster by its boundary. 

Thus, to deal with the distribution function of Eq. |T|), 
the relation between perimeter and area of clusters 
should be studied. It might be natural to suppose that 
the perimeter-area relation obeys the well known geo- 
metric formula: P cx vA However, it would be a very 
rough approximation, because this relationship holds for 
Euclidean geometric objects only. As was first found in 
Ref.— , the normal phase clusters can have fractal bound- 
aries, which exhibit non-Euclidean features. The fractal 
nature of such clusters exerts an appreciable effect on the 
dynamics of a magnetic flux in superconductorsi 13 i 40 i 41 



FIG. 3: The first six steps of the Cantor dust construction. 



IV. A BRIEF INTRODUCTION TO FRACTAL 
GEOMETRY 

The notion of a fractal as an object of fractional di- 
mension was first introduced by Mandelbrot 42 in 1967 
and since then it has received a lot of applications in var- 
ious domains of sciences . 43 ' 44 ' 45 ' 46 This concept is closely 
connected to ideas of scaling and self-similarity. Self- 
similarity is invariance with respect to scaling; in other 
words, an invariance relative to multiplicative changes 
of scale. Whereas a usual periodicity is invariance with 
respect to additive translations, self-similarity is a peri- 
odicity in a logarithmic scale. 

The simplest examples of self-similar objects are pre- 
sented by so called Cantor sets. Such a set has measure 
zero and, at the same time, has so many elements that 
it is uncountable. The Cantor set may be constructed 
in the following way. Let us draw a line segment, and 
erase its middle third. Then eliminate middle third of 
each remaining part, and so on. The first six steps of 
this procedure are shown in Fig. [3J The resulting set 
obtained by endless erasing of the middle thirds of re- 
maining intervals is called "Cantor dusf'4 3 - The Cantor 
dust forms a self-similar set: if one would take any line 
in Fig. [31 leave out the left or right half and magnify the 
remainder three times, this will result in the line segment 
immediately above it. So the Cantor dust is invariant to 
scaling by a factor s = 3. This set is not only a self- 
similar one, but it is a fractal as well. Let us look at 
that part of the set that falls within one s-th part of 
the original set, where s is a scaling factor, and ask what 
fraction of the set falls into that portion? In other words, 
how many subsets, each is similar to the original set, are 
there if the length is subdivided into s parts? The num- 
ber of such subsets is equal to N = 2, so one-half of 
the original set falls into one third of the initial length. 
The fractal dimension is defined as the logarithm of the 
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FIG. 4: Smoother Koch island: construction (a), and the final 
result (b). 



subset number divided by the logarithm of the scaling 
factor: D = bxN/bxs = In 2/ In 3 = 0.631. This formula 
represents the relation between the subset number and 
the scaling factor: N = s D . Thus the fractal dimension 
of the Cantor dust is less than unity. This fact reflects 
its dust-like consistency compared to a usual line. An 
Euclidean line is such a set that if we change the length 
scale, we recover the same set of points. Hence the fractal 
dimension of a line coincides with the topological dimen- 
sion, and both of them are equal to unity. Generally, a 
fractal object has a fractional dimension. A fractal set is 
such a set for which the fractal dimension strictly exceeds 
its topological dimension^ 

Another example of the fractal set is so called smoother 
Koch island^ 3 - Figure Ufa) shows the simple hexagon lat- 
tice formed by a hexagon surrounded by six identical 
hexagons. Let us note that each side of the big hexagon 
is broken up into three straight-line segments, each of 
length equal to \j\J~7 . The total area of seven hexagons 



is seven times the area of a little hexagon, whereas the 
perimeter of the entire figure is three times that of a little 
hexagon. Note that the total figure is not a hexagon so 
the figure is not similar to its parts. Although hexagons 
cover the whole plane, they do not form a hexagon again. 
But this peculiarity can be fixed up by making all the fig- 
ures to be similar. The first step is to make each straight 
line in Fig. [3Ja) to be similar to one third of the perime- 
ter of the entire figure. For this purpose we break up 
each side of a little hexagon into three segments again, 
each of the length of 1/7 of the side of the big hexagon. 
By repeating this breaking process endlessly, we will ar- 
rive at a figure that is sketched in Fig. BJb). The entire 
figure is similar to its seven parts. The obtained figure 
is a fractal called "smoother Koch island". Its scaling 
factor is equal to s = \fl by our construction, and the 
corresponding number of subsets is equal to TV = 3, so 
its fractal dimension is D = h\3/\xi^/7 = 1.129. The 
non-integer magnitude of the obtained value is consis- 
tent with the generalized Euclid theorem about measures 
of similar figures, which states that the ratios of corre- 
sponding measures are equal when reduced to the same 
dimension^ Therefore, it follows that 

P l ' D oc A 1 / 2 (11) 

which is valid both for Euclidean figures and for the frac- 
tal ones. 

Both the Cantor dust and the smoother Koch is- 
land are examples of determinate fractals: they are self- 
similar by deterministic construction. But the same ap- 
proach can be applied to the stochastic fractals, which 
are statistically self-similar. 43 The fundamental feature 
of any fractals, both determinate and the stochastic ones, 
is that its characteristic measures obey a scaling law 
that includes an exponent called the fractal dimension. 
In the following context, this is the scaling relation of 
Eq. (|lip where the perimeter and the enclosed area are in- 
volved, and the fractal dimension characterizes the clus- 
ter perimeter (so-called coastline dimension^ 3 .) 

V. FRACTAL GEOMETRY OF NORMAL 
PHASE CLUSTERS 

Using the relation of Eq. (jTTJ) between the fractal 
perimeter and the area of the cluster, as well as formula 
of Eq. (HHJ), we can get the following relation between 
the critical current of the cluster and its geometric size: 
I = aA~ D / 2 , where a is the form factor, and D is the 
fractal dimension of the cluster perimeter. In the gen- 
eral way the cluster area distribution can be described 
by gamma distribution ; 40 ' 41 which has the following cu- 
mulative probability function: 

W (A) = (r (g + I))" 1 7 (s + l^) (12) 

where T (is) is Euler gamma function, 7 (y) is the incom- 
plete gamma function, Aq and g are the parameters of 
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gamma distribution that control the mean area of the 
cluster A = (g + 1)Aq and its variance <j\ = (g + 1) Aq. 

The case of gamma distribution of the cluster areas 
is of most interest for optimizing the cluster structure of 
the composite superconductors, because this distribution 
has two independent parameters that can be varied in the 
cause of the film growth. One of the aim of such an opti- 
mization is to get the highest current-carrying capability 
of a superconductor. 47 ! 48 

In accordance with starting formulas of Eq. {]]) and 
Eq. gamma distribution of cluster areas of Eq. (Ti"2"|) 
gives rise to the critical current distribution of the form: 



F 



(z) = (r(.g + i))- 1 r(. 9 + i,«- 2 / 13 ) (13) 



where 



G = 



99+ 1 - (D/2)exp(0)r(. 9 + l, 



9 — g+ l + D/2, T(v, z) is the complementary incomplete 
gamma function, i = I / I c is the dimensionless electric 
current, I c = a(AoG) D ^ 2 is this the critical current of 
the transition into a resistive state. The found cumula- 
tive probability function of Eq. (flU)) allows us to derive 
the probability density / (i) = dF/di for the critical cur- 
rent distribution: 



2G^ 



DT(g + 1) 



-(2/D)(g+l) 



- 1 exp (-Gi- 2 ' D ) (14) 



This distribution allows us to fully describe the effect 
of the transport current on the trapped magnetic flux 
taking into account the fractal properties of the normal 
phase clusters. 

In order to clear up how the developed approach can 
be used in practice, the geometric probability analysis of 
electron photomicrographs of superconducting films was 
carried out. For this purpose electron photomicrographs 
of YBCO film prepared by magnetron sputtering have 
been scanned. The normal phase has occupied 20% of 
the total surface only, so the transport current can flow 
through the sufficiently dense percolation superconduct- 
ing cluster. The perimeters and areas of clusters have 
been measured by covering their digitized pictures with a 
square grid of spacing 60 x 60 nm 2 . The results of the sta- 
tistical treatment of these data are presented in Table U 
as well as in Fig. [5] The primary sampling has contained 
528 normal phase clusters located on the scanned region 
of a total area of 200 /im 2 . The distribution of the cluster 
areas is fitted well to exponential cumulative probability 
function, as is shown on the histogram in the inset of 
Fig. El The number of clusters that fall within the as- 
signed rank is plotted on the ordinate of this graph; the 
rank number is plotted on the abscissa. A high skewness 
(1.765) as well as the statistically insignificant (5%) dif- 
ference between the sample mean area of the cluster and 
the standard deviation also attests that there an expo- 
nential distribution of the cluster areas holds true. This 




0.001 



FIG. 5: Perimeter-area relationship for the normal phase clus- 
ters with fractal boundary. The solid line indicates the least 
square regression line, two dotted lines display the range of 
slope that the perimeter-area curves can have for any possi- 
ble fractal dimension D. (D = 1 - for clusters of Euclidean 
boundary, D = 2 - for clusters of boundary with the maxi- 
mum fractality). The inset represents the histogram of the 
cluster area sampling with rank spacing of 28680 nm 2 . 



distribution is a special case of gamma distribution for 
which g = 0, so the cumulative probability function of 
Eq. (fT2"|) can be simplified to the following form: 



W(A) = 1 - exp 



(15) 



The exponential distribution has only one characteristic 
parameter - the mean cluster area A. The obtained data 
allow us to find the perimeter-area relation for the nor- 
mal phase clusters. All the points in Fig. [5] fall on a 
straight line in double logarithmic scale with correlation 
coefficient of 0.929. Accordingly to the scaling relation 
of Eq. (jlip . the slope of the regression line gives the esti- 
mate of fractal dimension of the cluster perimeter, which 
is equal to D = 1.44 ± 0.02. The graph in Fig. O shows 
that the scaling law for perimeter and area, which is in- 
herent to fractals, is valid in the range of almost three or- 
ders of magnitude in cluster area. This relation between 
perimeter and area is consistent with the generalized Eu- 
clid theorem, so we can see that the ratios of perimeters 
and areas are equal when reduced to the same dimension. 
The scaling perimeter-area behavior means that there is 
no characteristic length scale in the range of two orders 
of linear size of the normal phase cluster. Whatever the 
shape and size of the clusters may be, all the points fall 
closely on the same straight line in logarithmic scale; so 
that there are no apparent kinks or bends on the graph. 

This point that the found value of the fractal dimen- 
sion differs appreciably from unity engages a great at- 
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Mean A, nm 2 


76540 


Sample standard deviation of A, nm 2 


72620 


Min value of A, nm 2 


2069 


Max value of A, nm 2 


401500 


Mean P,nm 


1293 


Sample standard deviation of P, nm 


962 


Alin value of P. nm 


Of; 


Max value of P, nm 


5791 


Correlation coefficient 


0.929 


Estimated fractal dimension D 


1.44 


Standard deviation of D 


0.02 



TABLE I: Statistics of normal phase clusters 



tention. What this means is that the fractal properties 
of the cluster boundary are of prime importance here. 
Two dotted lines in Fig. [5] bound the range of the slopes 
that the dependencies of the perimeter on the cluster 
area can have for any arbitrary fractal dimension. The 
least slope (upper line) corresponds to Euclidean clus- 
ters, which have the fractal dimension equal to topologi- 
cal one (D = 1), the most one (lower line) relates to the 
clusters of the greatest possible fractal dimension, which 
is equal to the topological dimension of a smooth surface 
(D = 2). Such a fractal dimension is inherent, for exam- 
ple, in Peano curves, which fill the whole plane.— What- 
ever the geometric morphological properties of clusters 
may be, the slope of their perimeter-area graphs will be 
always bounded by these two limiting lines. The found 
dependence (D = 1.44) runs just between them. 

The geometric probability properties of the normal 
phase clusters are responsible for main features of the 
critical current statistical distribution. Now, knowing the 
fractal dimension of the cluster boundaries, the change 
in the trapped magnetic flux caused by the transport 
current can be found with aid of Eq. @. The expo- 
nential distribution of the cluster areas of Eq. (fl"5j) gives 
rise to the exponential-hyperbolic distribution of critical 
currents 

^ )= exp(-p±^) 2/D+ Y 2 /^ (16) 

which follows from Eq. (flU)) at g = 0. Now the criti- 
cal current of the resistive transition I c , which appears 
in the expression for the dimensionless electric current 
i = 1 1 I c , can be found from a simpler formula: I c = 

(2/ (2 + £>))( 2+D )/ 2 a (Ay D/2 . The effect of a transport 
current on the trapped magnetic flux is illustrated in 
Fig. [6] for the case of Euclidean clusters (dotted curve) 
as well as for the clusters of found fractal dimension 
D = 1.44 (solid curve). 

In order to get the relationship between the dynamics 
of the trapped magnetic flux and the geometric morpho- 
logical properties of the superconducting structure the 
empirical function of the distribution of the critical cur- 
rents F* — F*(i) has been found. First, the empirical 
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FIG. 6: Effect of a transport current on the magnetic flux 
trapped in fractal clusters of a normal phase. The solid line 
shows the decrease in trapped flux for the fractal clusters of 
coastline dimension D = 1.44; the dotted line corresponds to 
the case of Euclidean clusters (D — 1); step line is the sample 
empirical function of critical current distribution. 

distribution function W* — W* (A) for the sampling of 
the areas of the normal phase clusters has been obtained. 
The value of W* (A) was calculated for each order statis- 
tic as the relative number of clusters of area smaller than 
a given value A. Next, the empirical distribution of the 
critical currents was computed for the same order statis- 
tics using the formulas: 



F* = 1 - W* 

(2±D\^+D)/2 (a\ 




This function, shown in Fig. [6] by the step line, gives a 
statistical estimate of the cumulative probability function 
of Eq. (|16p . As is seen from this figure, both distribu- 
tions coincide well in the range of currents i < 6. Start- 
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ing with the value of the current i = 6 the crossover from 
the fractal regime to the Euclidean one is observed. This 
transition into the Euclidean regime is over at large trans- 
port currents, when the magnetic flux changes mainly for 
the breaking of the vortices away from the small clusters 
(as the smaller clusters have the larger pinning force). 
The observed crossover has its origin in the finite reso- 
lution capability of measuring the cluster sizes. When 
estimating the fractal dimension, we have to take into 
account that the resolution of the measurement of any 
geometric sizes is finite. The peculiarity of the topolog- 
ically one-dimensional fractal curve is that its measured 
length depends on the measurement accuracy^ In our 
case such a fractal curve is represented by the boundary 
of the normal phase cluster. That is why just the statis- 
tical distribution of the cluster areas, rather than their 
perimeters, is fundamental for finding the critical current 
distribution. The topological dimension of perimeter is 
equal to unity and does not coincide with its fractal di- 
mension, which strictly exceeds the unity. Therefore the 
perimeter length of a fractal cluster is not well defined, 
because its value diverges as the yardstick size is reduced 
infinitely. On the other hand, the topological dimension 
of the cluster area is the same as the fractal one (both 
are equal to two). Thus, the area restricted by the fractal 
curve is a well-defined finite quantity. 

Taking into account the finite resolution effect, the 
perimeter-area relationship of Eq. (jlip can be re-written 
in the following form 



P(S)ocd 1 - D {A (5)) 



,D/2 



(17) 



where S is the yardstick size used to measure this length. 
This relation holds true when the yardstick length is 
small enough to measure accurately the boundary of all 
smallest clusters in sampling. When the resolution is de- 
ficient, the Euclidean part of the perimeter length will 
dominate the fractal one, so there is no way to find the 
fractal dimension using the scaling relation of Eq. (fl7|) . It 
means that if the length of a fractal curve was measured 
too roughly with the very large yardstick, its fractal prop- 
erties could not be detected, and therefore such a geomet- 
ric object would be manifested itself as Euclidean one. It 
is just the resolution deficiency occurs at the crossover 
point in Fig. [6] Starting with the cluster area less than 
some value (which corresponds to the currents of i > 6) 
it is impossible to measure all " skerries" and " fjords" on 
the cluster coastlines, whereas all the clusters of area less 
than the size of the measuring cell (3600 nm 2 that relates 
to the currents of i > 23), exhibit themselves as objects 
of Euclidean boundaries. This resolution deficiency can 
be also observed in Fig. [5j some points at its lower left 
corner are arranged discretely with the spacing equal to 
the limit of resolution (60 nm), because some marks for 
smallest clusters coincide for the finite resolution of the 
picture digitization procedure. 

The fractal dimension was found above by means of 
regression analysis of the whole primary sampling, where 
the very small clusters of sizes lying at the breaking point 



of the resolution limit were also included. In order to 
evaluate how the finite resolution affects the accuracy of 
the estimation of the fractal dimension, all the points, 
for which the resolution deficiency was observed, were 
eliminated from the primary sampling. So the truncated 
sampling has been formed in such a way that only 380 
clusters, for which the resolution deficiency is not ap- 
peared yet, have been selected from the primary sam- 
pling. The least squares treatment of these perimeter- 
area data gives the adjusted magnitude of the fractal 
dimension: D = 1.47 ± 0.03. The found value virtually 
does not differ from the previous one within the accuracy 
of the statistical estimation, whereas the correlation co- 
efficient (which becomes equal to 0.869) falls by six hun- 
dredth only. Therefore, we can conclude that the found 
estimate of the fractal dimension is robust. 

It is worthy of note that the above-described resolution 
deficiency refers only to the procedure of the cluster geo- 
metric sizes measurement. As may be seen from the Ta- 
ble U the characteristic sizes of the normal phase clusters 
far exceed both the coherence length and the penetration 
depth. Therefore, any effects of finite resolution related 
to the size of the vortices that scan the boundaries of the 
clusters in searching for weak links may be neglected. 



VI. FRACTAL SUPERCONDUCTING 
STRUCTURES IN A RESISTIVE STATE 

So, it has been revealed that the fractality of the cluster 
boundaries intensifies the pinning. As can be seen from 
Fig. [6] the decrease in the trapped magnetic flux at the 
same value of the transport current is less for larger frac- 
tal dimension. The pinning amplification can be charac- 
terized by the pinning gain factor 



fc # = 20 log 



A$ (D = 1) 
A$(L>) 



(18) 



which is equal to relative decrease (in decibels) in the 
fraction of magnetic flux broken away from fractal clus- 
ters of fractal dimension D compared to the case of Eu- 
clidean ones. The dependencies of the pinning gain on the 
transport current for different fractal dimension at g = 
are shown in Fig. [7] The highest amplification (about 
10 dB) is reached when the cluster boundaries have the 
greatest possible fractality. Figure [7] demonstrates that 
with increase in fractal dimension the trapped magnetic 
flux is changed less and less by the action of the trans- 
port current. The pinning gain of Eq. ()18[) character- 
izes the properties of a superconductor in the range of 
the transport currents corresponding to a resistive state 
(i > 1). At smaller current the total trapped flux re- 
mains unchanged (see Fig. [5]) for lack of pinning cen- 
ters of such small critical currents, so the breaking of 
the vortices away has not started yet. When the vor- 
tices start to leave the normal phase clusters and move 
through the weak links, their motion induces an electric 
field, which, in turn, creates the voltage drop across the 
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Eqs. (|T3|) , (|T6|l . so all the problems related to possible 
dependence of the flux flow resistance Rf on a transport 
current will not be taken up here. 

After substitution of the critical current distribution 
of Eq. (Ti"4|) into the expression of Eq. (fT9|), upon integra- 
tion, we get the final expression for the voltage-current 
( V-I) characteristics in the general case of the gamma- 
distribution of the cluster areas: 

£ = F<^T>(' r (» + ^- 2/D 

-G D / 2 T(g + l-^,Gi- 2 / D \\ (20) 



FIG. 7: Pinning gain for the different fractal dimension of 
the cluster perimeter in the case of exponential-hyperbolic 
distribution of the critical currents (g — 0) . 



sample. Therefore, the passage of electric current is ac- 
companied by the energy dissipation. As for any hard 
superconductor (that is to say, type-II, with pinning cen- 
ters) this dissipation does not mean the destruction of 
phase coherence yet. Some dissipation always accom- 
panies any motion of a magnetic flux that can happen 
in a hard superconductor even at low transport current. 
Therefore the critical current in such materials cannot 
be specified as the greatest non-dissipative current. The 
superconducting state collapses only when a growth of 
dissipation becomes avalanche-like as a result of thermo- 
magnetic instability. 

By virtue of the fact that any motion of the magnetic 
flux causes the energy dissipation in superconductors, the 
question of how such a process could be prevented, or 
only suppressed, is of prime practical importance. The 
study of resistive state peculiarities leads to the conclu- 
sion that the cluster fractality exerts influence on the 
electric held induced by the flux motion. The found dis- 
tribution of the critical currents allows us to hnd the elec- 
tric held arising from the magnetic flux motion after the 
vortices have been broken away from the pinning centers. 
Inasmuch as each normal phase cluster contributes to the 
total critical current distribution, the voltage across a su- 
perconductor can be represented as the response to the 
sum of effects made by the contribution from each clus- 
ter. Such a response can be expressed as a convolution 
integral 



V 
R~f 



(i-i')f(i')di' 



(19) 



where Rf is the flux flow resistance. The similar ap- 
proach is used universally in all the cases where the dis- 
tribution of the dcpinning currents takes place i 49 i 50 ' 51 i 52 
The subsequent consideration will be essentially concen- 
trated on the consequences of the fractal nature of the 
normal phase clusters specified by the distribution of 



In the special case of exponential cluster area distri- 
bution (g = 0) the general formula of Eq. (f2"0|) can be 
simplified: 



-C d ' 2 t{i-^Ci- 2 ' d ^ (21) 

where C= ((2 + D) /2) 2/D+1 . 

In extreme cases of Euclidean clusters and clusters of 
the highest fractality the expression of Eq. (I21|) can be 
further transformed: 

(i) Clusters of Euclidean boundary {D = 1 at g = 0): 



V 

S~ = t exp 



3^375 \ _ ^3 t 375 7rer f c 



where erfc (z) is the complementary error function, 
(ii) Clusters of boundary with the maximum fractality 
(D = 2 at g = 0): 



V 

— = i exp 
R f P \ , 



4Ei 



where Ei (z) is the exponential integral function. 

The V-I characteristics calculated using Eq. (f2"Tj) for 
several values of D are presented in Fig. [8] Two dot- 
ted lines bound the region the V-I characteristics can 
fall within for any possible values of fractal dimension. 
The solid curve represent the V-I characteristic of a su- 
perconductor with fractal clusters of previously obtained 
fractal dimension D = 1.44. The inset in Fig. [5] shows 
the region of resistive transition under magnification. As 
may be seen from this graph, the critical current i c is 
preceded by some initial region of the finite voltage drop 
starting with i on , so the resistive transition is not ab- 
solutely abrupt. The existence of this initial section on 
the V-I characteristic arises from the peculiarities of the 
fractal distribution in the range of small currents. It 
has been just a similar initial region of fractal dissipa- 
tion has been observed in high-resolution measurements 
of dynamical resistance of HTS- normal metal composites 
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FIG. 8: Voltage-current characteristics of superconductors 
with fractal clusters of a normal phase for exponential- 
hyperbolic distribution of the critical currents (g — 0). Dot- 
ted lines correspond to extreme cases of Euclidean clusters 
(D = 1) and clusters of the most fractality (D = 2); solid 
line - to superconductor containing the normal phase clusters 
of fractal dimension D = 1.44. The inset shows the region 
near the resistive transition in an enlarged scale. The initial 
dissipative range between i on and i c is clearly seen. 



(BPSCCO-Ag) as well as in polycrystalline YBCO and 
GdBCO samples.ii^ 

All the curves of V-I characteristics are virtually start- 
ing with the transport current value equal to unity. When 
the current increases the trapped flux remains unchanged 
until the vortices start to break away from the pinning 
centers. As long as the magnetic flux does not move, no 
electric field is arisen. Figure [5] shows that the fractal- 
ity reduces appreciably an electric field arising from the 
magnetic flux motion. This effect is especially strong in 
this range of the currents (1 < i < 3), where the pin- 
ning enhancement also has a maximum (see also Fig. [7|) . 
Both these effects have the same nature, inasmuch as 
their reason lies in the peculiarities of the critical cur- 
rent distribution. The influence of the fractal dimension 
of the cluster boundary on the critical current distribu- 
tion is demonstrated in Fig. [5J As may be seen from this 
graph, the bell-shaped curve of the distribution broad- 
ens out, moving towards greater magnitudes of current 
as the fractal dimension increases. It means that more 
and more of the small clusters, which can best trap the 
magnetic flux, are being involved in the game. Hence 
the number of vortices broken away from pinning centers 
by the Lorentz force is reducing, so the smaller part of a 
magnetic flux can flow. The smaller part of a magnetic 
flux can flow, the smaller electric field is created. In turn, 
the smaller the electric field is, the smaller is the energy 
dissipated when the transport current passes through the 
sample. Therefore, the decrease in heat-evolution, which 
could cause transition of a superconductor into a nor- 
mal state, means that the current-carrying capability of 
the superconductor containing such fractal clusters is en- 
hanced. 



FIG. 9: Influence of the fractal dimension of the cluster 
perimeter on the exponential-hyperbolic distribution of the 
critical currents (g = 0). The inset shows the dependence of 
the mean critical current i on the fractal dimension. 



The re-distribution of probability with increasing in 
the fractal dimension can be characterized by the statis- 
tical moments of the critical current distribution, namely, 
the mean critical current 



D/2 T(g + l-D/2) 

r(s + i) 



and the variance of critical currents 



G 



D 



T(g+1-D) {T(g+l 



r G? + i) 



1) J 



As may be seen from the inset in Fig. [9l the mean criti- 
cal current grows nonlinearly with increase in the fractal 
dimension that reflects the broadening of the "tail" of 
the critical current distribution. This "tail" is so elon- 
gated that the variance of critical currents becomes in- 
finite in the range of fractal dimensions D > g + 1. At 
the same time, the mode of the critical current distribu- 
tion is a finite function, which increases monotonically 
in the range of all possible values of fractal dimension: 
modef (i) — (G/8) D ^ 2 . The distributions with divergent 
variance are known in probability theory - the classic ex- 
ample of that kind is Cauchy distribution^ However, 
such an anomalous feature of the critical current distri- 
bution of Eq. is of special interest, inasmuch as the 
current-carrying capability of a superconductor would be 
expected to increase just in the region of giant dispersion. 
Then the statistical distribution of critical currents has a 
very broad "tail" containing the contributions from the 
clusters of the highest depinning currents. 



VII. CONCLUSION 

So, in the present work the fractal nature of the normal 
phase clusters is revealed, and relation between the frac- 
tal properties of the clusters and dynamics of the trapped 



12 



magnetic flux is established. The fractal distribution of 
the critical current is obtained. It is found that the frac- 
tality of cluster boundary strengthens the flux pinning 
and thereby hinders the destruction of superconductivity 
by the transport current. V-I characteristics of fractal su- 
perconducting structures in a resistive state are obtained. 
It is revealed that the fractality of the boundaries of the 
normal phase clusters reduces the electric field arising 
from magnetic flux motion, and thereby raises the criti- 
cal current of a superconductor. 

At the same time, the following problems are open for 
further development in this field. First, a central point of 
the presented approach is that the magnetic flux is held 
on the boundaries of clusters. In experiments with poly- 
crystalline magnesium diboride^ it has been found that 
in this new superconductor the similar situation takes 
place, but instead of the cluster boundaries, the grain 
boundaries act in the same way. Second, the fractal dis- 
tribution of the critical currents have anomalous statis- 
tical properties caused by a divergence of some its sta- 
tistical moments so this point is worthy for further in- 



vestigations. Then, the special cases of the distribution 
of entry points into weak links, which can be realized in 
anisotropic composites, are of interest. First of all, we 
mean the entry point distributions that take place in su- 
perconducting tapes and wires. Last, obtained results 
open the possibilities for increasing the critical current 
value of percolative superconductors by optimizing their 
geometric morphological properties, and enable to im- 
prove the technology of preparation of superconducting 
films of high critical current value. The study of these 
questions will allow us to get a better insight into the vor- 
tex dynamics in the superconductors containing fractal 
clusters. 
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